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Abstract 



Magnetic Reconnection is an efficient and fast acceler- 
ation mechanism by means of direct electric field accel- 
eration parallel to the magnetic field. Thus, accelera- 
tion of particles in reconnection regions is a very im- 
portant topic in plasma astrophysics. This paper shows 
that the conventional analytical models and numerical 
test particle investigations can be misleading concern- 
ing the energy distribution of the accelerated particles, 
since they oversimplify the electric field structure by the 
assumption that the field is homogeneous. These inves- 
tigations of the acceleration of charged test particles are 
extended by considering three-dimensional field config- 
urations characterized by localized field-aligned electric 
fields. Moreover, effects of radiative losses are discussed. 
The comparison between homogeneous and inhomoge- 
neous electric field acceleration in reconnection regions 
shows dramatic differences concerning both, the maxi- 
mum particle energy and the form of the energy distri- 
bution. 

PACS numbers: 52.20 ; 52.65.C ; D95.30 ; 95.30Q 



I. Introduction 



The acceleration of charged particles in space and as- 
trophysical plasmas is a matter of debate iii contem- 
porary plasma astrophysics (cf. reviewsH andcl). Beside 
tb^" standard" acceleration mechanisms as shock 



, Fermi acceleration!^, and magnetic turbulence (13,111) 
acceleration by magnetic field-aligned electric fields re- 
ceives more and more attention as it offers the capabil- 
ity for, efficient and fast energization of relativistic parti- 
cles (EI,u). The formation and amplification of such elec- 
tric fields in the context of magnetic reconnection pro- 
cesses has been discussed. ^E-difforent cosmical plasmas 
and parameter regimes (ll^l,EI,tII,ll3jl3). Since magnetic re- 
connection can be regarded as an unavoidable dissipation 



and relaxation process in magnetized plasmas under the 
influence of externally driven shear motions, the princi- 
pal importance of magnetic reconnection electric fields 
is now widely accepted. 

Analytical and numerical studies have signifi- 
cantly improved our understanding of energetic par- 
ticle J-raiectodes, an reconnectionpjfaafigurations with- 
out (yjil,y,El,Bjia) and with a magnetic field 

component parallel to the reconnection current sheet. 
Nevertheless, the particle trajectories and, in particu- 
lar, the actual energy gains of accelerated charged par- 
ticles in the reconnection current sheets prove to be still 
a crucial matter of debate, since the studies mentioned 
above have been performed for idealized analytical elec- 
tric and magnetic field configurations which differ signif- 
icantly from the three-dimensional configurations char- 
acterized by localized field-aligned electric fields as ob- 
served, e.g. in the jSolar corona (cf£3) and the Earth's 
rpagnetosphere (cf.cS). Test particle simulations of ion 
(Ej) and electron (c3) injection in a high-beta magne- 
tospheric substorm configurations have shown that the 
consideration of more complicated 3d field configurations 
rather than idealized 2D ones is crucial for the under- 
standing of particle acceleration in reconnection regions. 
In this contribution we study the acceleration of charged 
particles up to relativistic energies in reconnection re- 
gions. The electromagnetic fields of the reconnection re- 
gion are modeled by 3D-magnetohydrodynamic (MHD) 
simulations of a reconnection zone in a current carry- 
ing magnetic flux tube. Since the astrophysical evidence 
for the existence of relativistic particles relies only on 
the electromagnetic radiation emitted by of these parti- 
cles, we calculate the motion of high-energy particles by 
means of relativistic test particle simulations including 
radiative forces due to inverse Compton scattering and 
synchrotron radiation. We consider our study as a step 
forward to a self-consistent description of particle accel- 
eration in quite realistic reconnection configurations. 

The manuscript is organized as follows: The next 
section presents the motivation for our new approach. 
Sec. II] describes the electromagnetic field configuration 
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we used and how they are obtained from MHD simula- 
tions. These fields represent the external force field, to 
be supplemented by radiative forces, for the relativistic 
test particle calculations, described in Sec. IV. The re- 



sults of these test particle simulations are compared with 
the analytical and numerical X-type configuration of the 
reconnection magnetic field in Sec. Sec. VI is devoted 
to a discussion of our results. 



II. Motivation for our new Approach 

Particle acceleration in magnetic reconnection zones has 
been extensively discussed and studied analytifiaUv-aiid 
numerically (we refer to the analytical paperaijj23,E3,E3, 
and refexfiiaces therein); numerical treatments can be 
found inE£l,t3, and references therein). These investiga- 
tions rely on one common assumption The homogeneity 
of the electric field. The authors have used an analyt- 
ical formulation of the electric field configuration and 
therefore the calculations are highly idealized. The field 
configurations found in the papers mentioned above are 
nearly identical and can be written as 
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E{r)=Ea 



The process of reconnection in two dimensions is still 
subject of intense scientific investigation up to-, the 
present day. Besides first analytical approaches (pl;EJ) 
investigations of the much more complex 3D reconnec- 
tion process are restricted to numerical modeling. In a 
2D description the X-point configuration 



y 



B{r) = B, I f-^ 




(2) 



that describes a hyperbolic magnetic field with a neu- 
tral line in the center (B(0,0, z) = 0), has become the 
standard model of reconnection. With Ampre's Law 



4vr . „ „ 

— J = V X B, 

c 

Ohm's Law 

Tjj = E + v/cy^ B, 



(3) 



(4) 



and invariance in one dimension, say dz- = 0, one can 
easily derive the electric field 
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v/c X B. 



(5) 



If one assumes a constant resistivity and that the con- 
vective term v/c x B \s small as compared to the first 
term in (Equ. 0) in the localized reconnection zone the 
constant, homogeneous, infinitely extended electric field 
is given by 
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(6) 



Such an electric field is, without question, unphysical. In 
case of a 2D reconnection calculation this is not a severe 
problem, because the electric field is never used as a real 
physical quantity. For this reason one can argue that 
it is merely a mathematical auxiliary variable and only 
the magnetic field is needed to describe the complete 
electromagnetic configuration. 

However, this X-type field is a valid approximation of 
the inner part of a reconnection region around the neu- 
tral line. Only the application of this idealized config- 
uration in the context of particle acceleration gives rise 
to serious problems. The electric field is the source of 
energy for the accelerated particles and exactly this very 
reason for acceleration is, as stated above, unphysical. 
At this point, we should ask how sensible are the results 
concerning particle energy spectra, duration of acceler- 
ation etc., derived from an idealized, homogeneous elec- 
trical field? For example, there are particle trajectories 
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for which particles will gain infinite energy. This comes 
from the fact that the electric field lacks any dependence 
on the invariant direction. 

In order to accelerate p^articles efficiently a strong 
guiding component Bz (E3) is needed. With such a 
strong Bz component the magnetic field is almost paral- 
lel to the electric field in a region + By around 
the neutral line. In addition, this guiding component 
captures the particles in the zone, at least as long as 
their gyroradius r^c is smaller than the width of the re- 
connection zone. These two facts guarantee the efficient 
acceleration described incJ. On the other hand this con- 
figuration in the center of a X field with strong Bz resem- 
bles remarkably that of a simple capacitor. Any particle 
following the field lines is accelerated exactly linear, like 
in a capacitor. The energy gained by such a particle is 
directly given by 



AS = qEzAL. 



(8) 



A£ denotes the energy gain and q is the charge of the 
particle that is accelerated along the magnetic field line 
over the distance AL. 
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Fig. 1. Sketch of the generic initial configuration. 

The energy gain does not depend on the detail of the 
reconnection process, but on the geometry of the accel- 
eration region only. 

Moreover, one has to take into account the convec- 
tive electric field components and Ey, which are 
completely neglected by most authors, though they 
naturally arise in realistic descriptions of reconnection 
zones. Such fields are generated by external shear 
flows which are the macroscopic generator of electro- 
magnetic energy converted partly into particle kinetic 
energy in the localized dissipative reconnection zones. 
Calculations with three-dimensional mamptohydrody- 
namic (3D-MHD) codes clearly show (e.g.El) that those 
components are often much stronger than the fleld line 
parallel electric field . The precise influence of those 
components is not known, nevertheless we should not 
neglect them strait from the beginning, without know- 
ing anything about the strength of their influence. We 
speculate that such convective electric fields can filter 
out particles from an efficient acceleration due to their 
strength and structure and dependent on the initial po- 
sitions of the particles. In other words, the form of the 
energy distribution may be closely connected to these 
electric fields. 

The above questions convinced us to approach again 
this well known and well studied phenomenon of par- 
ticle acceleration in magnetic reconnection zones. Con- 
trary to other investigations, we use a field configuration, 
which is produced with the help of a full 3D-MHD code. 



III. The 3D-MHD Simulations of Magnetic Recon- 
nection 

The MHD framework for the particle simulations of 
electron acceleration in reconnection regions is— built 
up by means of a three-dimensional MHD codeE3. It 
is motivated by previous works on reconnection in 
the coBtestpOfrftlaHetary magnetosphere-ionosphere cou- 
pling (EilJljj23j23,o) and of active galactic nucleill3. How- 
ever, the situation we have in mind (cf. Fig. |l|) is quite 
general and may be regarded as a rather generic recon- 



nection configuration that arises when MHD shear flows 
working as electromagnetic generators sheared magnetic 
fields. For convenience, we start from a two-dimensional 
force-free magnetic field 



B = Byo tanh(a::) Gy — W 
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yo 



cosh (x) 



(9) 



We choose for the constant main component and the 
shear component of the magnetic field B^a — 5 and 
ByQ = 1 in normalized units. The plasma is chosen to 
be isothermal as well as homogeneous, initially. We note 
that, alternatively, an appropriate equilibrium that in- 
cludes thermal pressure forces could be realized. This 
ideal equilibrium is perturbed by 



tanh(2a;) 

Vy(x, z) = Vyo -j-T^ exp 

cosh (f) 



As a boundary condition for t > we use 



Vy{x,Z 



l) = VyQ 



tanh(2j:) 
cosh^ 



(10) 



(11) 



These velocity perturbation with an amplitude chosen as 
Vyo — 0.5% of the Alfven velocity mimics some external 
convective plasma motion which leads to a further shear 
of the magnetic field and an increase of the field-aligned 
current density. The associated magnetic perturbation is 
transported along the ^-component of magnetic field via 
shear Alfven waves. Here, we follow the general concept 
that plasma systems that are compelled by such external 
forces when some critical parameter(s) is (are) exceeded 
non-linearly develop into a state of lower energy thereby 
releasing stored free energy in form of dissipative plasma 
heating and particle acceleration. 




Fig. 2. Some characteristic magnetic fieldlines af- 
ter t = 120rA 
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Fig. 3. Snapshot of the x-, y- and z-component (from top 
to bottom) of the reconnection magnetic field at the height 
of the central reconnection region after t = 120rA. 

In the considered configuration the magnetic field- 
aligned current density plays the role of the critical pa- 
rameter. When the current density exceeds a critical 
value (we choose a marginal current density at the begin- 
ning) a current dependent resistivity, which is localized 
in the y- and z-direction is switched on. This violation 
of ideal Dhm's law starts the magnetic reconnection pro- 
cess (cf.El). Fortunately, for our purpose we have not to 
dwell on the specific |iiature of the finite resistivity, may: 
it be microturbulentEj dissipation or particle inertiaEB 
or whatever. For the MHD simulations we choose the 
following technical details: 
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Fig. 4. Snapshot of the x-, y- and z-component (from 
top to bottom) of the reconnection magnetic field at the 
y = 0-half-plane after t = 120ta. 

The dimensions of the numerical box are given by 
X G [-10,10], y e [0,30], and z G [0,60] in normal- 
ized units. We make use of line symmetry as a bound- 
ary conditions dX y — Q, i.e. the x- and y-components 
of the magnetic field and the velocities are chosen anti- 
symmetric whereas the z-components of these quantities 
are chosen symmetric about y = 0. Boundary condi- 
tions in X are chosen as symmetric and in z as well as 
at 2/ = j/max we use extrapolation. The simulations are 
carried out with 49 grid points in the x- direction, 39 
grid points in the y-direction and 105 grid point in the 
z-direction. 
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Fig. 5. Snapshot of the x-,y- and z-component (from top 
to bottom) of the reconnection electric field at the height of 
the central reconnection region after t = 120rA- 

The numerical grid is a non-uniform one with a max- 
imum resolution of 0.05 in the x-direction, 0.4 in the 
y-direction, and 0.2 in the z-direction. 

Some characteristic fieldlines are shown in Fig. || to 
give an idea of the overall structure of the magnetic field. 
In Figs. ^ to ^ we show the Cartesian components of the 
electric and magnetic fields which are chosen as an in- 
put for the particle simulation studies. The figures show 
the fields after t = 120ta where ta is the characteris- 
tic Alfven velocity of the system. At this state of the 
temporal evolution the reconnection process can be con- 
sidered as quasi-stationary, i.e. the magnetic flux is as 
rapidly piled up by the shear flow as it is annihilated by 
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Fig. 6. Snapshot of the x-,y- and 2-component (from 
top to bottom) of the reconnection electric field at the 
y = 0-half-plane after t = 120ta. 

magnetic reconnection. 

The x-component of the magnetic field results from 
the shear fiow and from the divergent reconnection flow 
(Fig. ^, top panel) and (Fig. ^, top panel). The y- 
component of the magnetic field (middle panels of Fig. ^ 
and Fig. ^ is associated with the current sheet and mod- 
ified by reconnection convergent flows {vx{x)). The z- 
component (lower panels of Fig. ^ and Fig. ^) shows a 
rather moderate variation due to the reconnection pro- 
cess. The three-dimensional electric vector field at the 
height of the central reconnection region is illustrated in 
Fig. ^ The convective x-component (top panels of Fig. ^ 
and Fig. ^ is the result of both the applied sheared flow 
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and the divergent reconnection flow whereas the convec- 
tive y-component (middle panels of Fig. || and Fig. ||) is 
due to the convergent reconnection flow. The convective 
electric fields are about one order of magnitude higher 
than the ^-component of the electric field (lower pan- 
els of Fig. H and Fig. ||) which is caused by the finite 
anomalous electric resistivity Ez ^ i]jz ~ E^^ w rjj^^ . E^ 
is localized in the x- and y-direction (we have chosen for 
the resistivity profile 77 ~ cosh{y{iv)/2.)) and is most 
important for the acceleration of charged particles. 

It is mainly responsible for the non-vanishing pJp*i- 
gated generalized field-aligned electric potential (□,EII) 
U = — J E^^ds (where the integration is extended over 
the arch length s). Such potentials are observed, e.g. as 
"aurorjal potential structures" in the Earth's magneto- 
sphereEa which are responsible for the acceleration of au- 
roral electrons. The three-dimensional electric and mag- 
netic vector field are used as an input for the test particle 
simulation studies of relativistic electrons as discussed 
in the next chapter. Note, that these rather complex 
fields, which are due to shear Alfven waves, shear mo- 
tion and reconnection dynamics, differ significantly from 
idealized configurations usually discussed in the litera- 
ture. We consider the electromagnetic field shown as a 
quite typical reconnection configuration in systems that 
are characterized by shear flows. Since it represents the 
self-consistent reconnection dynamics it can be qualified 
as a realistic configuration in comparison to idealized 
two-dimensional ones. 



IV. Particle Acceleration 



We use the fields described in Sec. Ill for electron ac- 
celeration. By means of a fully relativistic test particle 
code, including energy losses due to synchrotron emis- 
sion (SY) and inverse Compton scattering (IC) we ana- 
lyze the fate of the relativistic particles. The term "test 
particle" in this context means that any back-reaction of 
particles on the electromagnetic field configuration is ne- 
glected. As long as only a fraction of particles is actually 
accelerated, i.e. as long as the particles leave the acceler- 
ation zone on time scales much shorter than the Alfvenic 
time this is a very good approximation, which is always 
fulfilled during our simulations, since the particles are 
efficiently accelerated along the magnetic field line and 
move with almost the velocity of light. Moreover, no 
interaction between the particles is included. 

For the integration of the relativistic equations of mo- 
tion 



dr dp 

— ^v -i- = F 

dt dt 
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(12) 



we use the classical fourth order Runge-Kutta algorithm 
together with an adaptive stepsize controlo. Our 3D- 
MHD fields are only known on a discrete 3D grid. For 
this reason, we have to interpolate these fields for posi- 
tions in between by means of a simple linear interpola- 
tion algorithm. 

As described in Sec. 



Ill the electromagnetic fields 



can be regarded as quasi-stationary after some hundred 
Alfvnic times t ^ 120rA. Moreover, the test electrons 
are-accelerated on a very short time scale tacc — 10~^ta 
(cf.^). For this reason there is no need to consider the 
precise dynamic evolution of the MHD fields. We can use 
them as static background fields without any temporal 
changes and this is also a very good approximation. In 
addition to the Lorentz force i^L = q (|^+ "^/c x B) our 
code considers the complete loss forceEJ 
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(13) 



where q, m, v and 7 denote the charge, the mass, the ve- 
locity and the Lorentz factor of a particle, respectively. 
This loss force describes all particle-field interactions, in- 
clusive synchrotron radiation (SY) and inverse Compton 
scattering (IC). SY presents the radiation of relativis- 
tic electrons, gyrating around the magnetic field, i.e. it 
influences the particle momentum perpendicular to the 
magnetic field. If an acceleration process injects parti- 
cles merely along the field lines these particles undergo 
no synchrotron losses. However, since we inject parti- 
cles isotropically from all directions, the particles will 
have flnite pitch angles at the beginning. Even for very 
anisotropic energy distributions in the sense that the 
particles move only parallel to the magnetic fleld lines 
particles may suffer radiative losses caused by inverse 
Compton scattering, if an intense photon bath is present 
in the environment of the acceleration zone. The radia- 
tion mechanisms depend on the energy of the particles 
and the energy density of the magnetic field Uung oc B'^ 
(SY) or the photons C/Rad (IC), respectively. 

Qualitatively, it is clear that relatively strong (in terms 
of efficiency and short time scales) particle acceleration 
via magnetic reconnection requires significant magnetic 
field strengths. Since astrophysical plasmas are highly 
conducting media, strong magnetic fields are related 
with either fast plasma motions or high plasma pressure 
and temperature, i.e intense radiation fields. Such condi- 
tions are typical, e.g., for compact astrophysical objects, 
as neutron stars or black holes and magnetospheres of 
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young stellar objects. Often the magnetic field energy 
density and the radiation energy density are of compa- 
rable order, or at least proportional to each other. This 
quality is again connected with the high electric conduc- 
tivities of astrophysical plasmas on macroscopic spatial 
scales. If plasma is accelerated due to some radiation 
pressure, the magnetic field is strongly coupled to the 
plasma motion and is amplified by field line compression 
and twisting. The plasma temperature is directly related 
to the radiation temperature. Thus, a significant varia- 
tion of the emitted radiation will consequently produce 
varying plasma temperatures and pressures ,|jyhich fi- 
nally leads to magnetic field amplification (e.g.c3). Thus, 
the energy losses of the accelerated particles can be han- 
dled by a simple scaling Uuag t/rad- 



V. Comparison 

We use our code for the calculation of electron energy 
spectra which one would expect to be the final energy 
distributions in some selected configuration. In order 
to perform this we calculate the trajectories of more 
than 8000 electrons, which start at random places with 
a random momentum and determine their final energies, 
when they leave the numerical bounding box of the MHD 
fields. The initial positions are equally distributed in the 
whole reconnection region. The initial momentum of the 
electrons is chosen isotropic with absolute values that 
correspond to a Maxwellian distribution of the initial 
energy. 
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Fig. 8. Final energy distribution, calculated with realis- 
tic 3-D-MHD fields. 

In Fig. 1^ the initial energy distribution of the elec- 
trons is shown. Here and in all following plots "Rela- 
tive Frequency" is the normed (J dE = 1) percentage 
of particles per energy interval. All runs start from this 
distribution. Fig. ^ shows the final energy distribution 
of electrons accelerated in our quite realistic 3D-MIID 
fields (see Sec. III). Most electrons are accelerated due 
to the strong guiding component Bz- Energies in the 



range from 10 



to nearly 10 moC are seen. One 



first remarkable fact is the existence of electrons nearly 
without any energy gain. A second interesting feature is 
the obvious depression in the middle of this spectrum. 
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Fig. 7. Initial energy distribution (common to all simu- 
lations) . 
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Fig. 9. Final energy distribution, calculated with re- 
alistic 3-D-MHD fields and consideration of loss processes. 

f/Rad ~ 10 • (7Mag. 
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Fig. 10. Final energy distribution, calculated with a sim- 
plified, homogeneous electric field and the realistic magnetic 
field. 

It seems as if there exists some sort of selection pro- 
cess, which does not allow a particle to achieve an energy 
between IO^ttioC^ and \Q'^mc<? . This might indicate that 
there are unstable trajectories, which lead exactly to this 
energy range. The high energetic end of the particle dis- 
tribution is determined by the thickness of the resistive 
region. Only in the region of finite resistivity a parallel 
electric field is present, which can accelerate the particles 
efficiently. 
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Fig. 12. Pitch angle distribution of the electrons simu- 
lated in the oversimplified electric field. 

The gyro radius of an accelerated electron is given by 

VTcimeC 



eB 



(14) 



which depends on the energy of the electrons (uto is the 
thermal velocity of the electrons). Particles are accel- 
erated, in the best case, as long as their gyro radii are 
smaller than the extent of the resistive region. This cri- 
terion gives the maximum energy gain achievable. The 
effect of considering the loss force (Equ. |l^) is shown in 
Fig. 1^. In order to illustrate those effects we overestimate 
the strength of the loss force by assuming an extremely 
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Fig. 11. Pitch angle distribution of the electrons simu- 
lated in the realistic 3D-MHD fields. 
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intense photon bath with an energy density ?7Rad of 
about ten times the maximum magnetic energy den- 
sity Uyiag- Very few electrons are now able to achieve 
more than W^rrieC^. This is the energy, where the loss 
of energy due to radiative processes balances the energy 
gain from the electric field. The characteristic features 
mentioned before remain unchanged apart from the new 
upper limit for the final energy, which leads to a "pile- 
up" at about lO^mcC^. 

After all, it can be stated that the usage of rather re- 
alistic 3D-MHD fields, i.e. a field configuration that rep- 
resent evolved reconnection electromagnetic fields rather 
than analytical idealized fields, gives a rich spectrum dis- 
playing a great variety of features. This changes drasti- 
cally in the case of an oversimplified homogeneous elec- 
tric field. The next example is calculated with the elec- 
tric field 



E{r) 



(15) 



where Ez^max denotes the maximum of the z component 
of the numerical electric field used before. The resulting 
energy distribution, calculated in the same way as be- 
fore, is shown in Fig. |l0|. We should point out that for 
this simulation nothing has been changed in the mag- 
netic field configuration. Fig. 00^ shows a very narrow 
spectrum around the energy KTVicC^ completely with- 
out any feature so characteristic for the spectrum of 
Fig. ^. First of all, we can clearly see that every elec- 
tron is accelerated. This unphysical result comes at no 
surprise, because a constant electric field E^ in combi- 
nation with the strong guiding component implies 
an infinite resistive region. The obvious absence of any 
structure shows, that the origin of the depression has 
to lay mainly in the electric field. An oversimplification 
of E neglects this interesting effect. The enlargement 
of the resistive region to the whole numerical box also 
erases the upper limit of acceleration. The maximum 
energy is given by the potential difference between on 
side of the box and the other, but that is a pure trivial 
geometric limit. 

In addition, there is no similarity between the rather 
realistic spectrum Fig. |^ and that calculated with a ho- 
mogeneous electric field Fig. |lo[ 

Comparison of the pitch angle distribution (Fig. ^l]) of 
the realistic simulation with the one of the oversimplified 
electric field simulation (Fig. |l^) also shows a striking 
difference. While both distributions clearly show that 
a vast majority of electrons is moving mainly along the 
magnetic field lines, again a broad spectrum over tree 
decades and a small amount of electrons, which have a 
pitch angle greater than ^ evolves (Fig. O). The pitch 



angle distribution in Fig. on the other side is lacking 
any hint of the electron population with an pitch angle 
greater than « 0.01. An homogeneous electric field can, 
in fact, not be regarded as a valid approximation of a 
realistic field configuration. 

At this point the question arises, how important is the 
structure of the magnetic field for the final spectra any- 
way? In order to answer this question we finally consider 
the spectrum for 
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This spectrum (Fig. |T^) shows a striking similarity 
with the previously shown one Fig. |l^. Both spectra are 
dominated by the geometry of the numerical box. All 
of the features seen in Fig. ^ are coming from the com- 
plicated structure of the electric field. It seems that the 
structure of the magnetic field serves only as " stage" for 
the electrons without significant influence on the spectra. 
Nevertheless, we should keep in mind, that the B field 
determines the structure of the E field. 



VI. Discussion 

Astrophysical plasmas like the interstellar or intergalac- 
tic medium and compact systems like accreting stellar 
objects are magnetized and highly conducting collision- 
less media under the influence of external forces. The 
main driving forces for plasma-magnetic field motions 
are gravity, rotation and pressure gradients. Due to the 
lack of significant particle collisions astrophysical plas- 
mas are easily forced to steep magnetic field gradients, 
which ultimately will evolve into relaxed field configura- 
tions via magnetic reconnection. In that sense magnetic 
reconnection is an unavoidable result of external distor- 
tions in coUisionless magnetized plasmas. The important 
observable consequences of reconnection for astrophysi- 
cal applications follows from the fact that it accounts 
for heating and particle acceleration by locally convert- 
ing the stored magnetic energy into bulk heating of the 
plasma via Ohmic heating or acceleration of a popula- 
tion of high energy particles via magnetic field-aligned 
electric fields that are associated with the finite conduc- 
tivity. 

The heating capabilities of magnetic reconnection 
have been investigated in detail, since the heating rate 
(X rip is a macroscopic quantity which depends only on 
the density and temperature, but not on the details of 
the single particle response to the electromagnetic con- 
figuration inside a reconnection zone. In that sense it 
behaves like a macroscopic transport property. 
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The acceleration of particles instead depends on the 
local properties of the reconnection region. The ener- 
gization of a charged particle above all depends on the 
effective electric force the particle "feels". Depending 
on its initial position and initial energy the particle en- 
ergy gain differs. Since the only source of information 
of high energy astrophysics is the electromagnetic radi- 
ation, the energy distribution of relativistic particles is 
one of the most required issues concerning the interpre- 
tation of nonthermal radiation sources. Thus, for any 
acceleration process discussed it is most important to 
investigate the resulting energy distribution functions in 
dependence on the various parameters on which the ac- 
celeration process depends on. In this context, the two- 
dimensional models published so far for the acceleration 
of charged particles by magnetic reconnection are not 
satisfying, since they describe acceleration in linear con- 
ductors. They do not describe the effect of the finite 
extension of the acceleration zone and show always a 
purely monoenergetic distribution at the maximum en- 
ergy the particles are "allowed" to reach. The maximum 
energy in these models is simply given by the product of 
the homogeneous electric field times the chosen length 
of the reconnection zone, which is supposed to be a free 
parameter. This is definitely incorrect; the system size 
must be finite and must depend on the process respon- 
sible for the violation of ideal conductivity. 

We have shown the pitfalls of the two-dimensional 
ansatz and present 3D-simulations of MHD-rcconncction 
plus test particle simulations within the MHD-frame, in- 
cluding radiation losses (synchrotron and inverse Comp- 
ton scattering) . The precise structure of the electric field 
is the determining factor of the acceleration process in 
the MHD reconnection region (cf. Sec. 0). One has to be 
extremely careful in simplifying the electric field; in our 
example the spectrum calculated with a homogeneous 
electric field is no good approximation to the one that 
results from a physically motivated field. 

Since the investigated magnetic field configuration is 
a valid 3D generalization of a 2D X-type field the us- 
age of an homogeneous electric fields in X-type electric 
fields should also be misleading as far as the spectrum is 
concerned. 

Whereas particle acceleration in a simplified X-type 
configuration with a homogeneous electric field is an in- 
teresting subject of it's own, it is obviously not appro- 
priate for the study of particle acceleration in magnetic 
reconnection zones. We have shown that for this task 
we need much more realistic fields and at least numeri- 
cal approaches in future should use 3D-MHD fields. A 
further extension of the work presented in this contri- 
bution will be the consideration of multiple XJine and 
0-point reconnection configurations (c.fo andcj for the 
2D case) which should show a larger class of different 



kinds of particles populations. 

The strong infiuence of radiative losses on the resulting 
energy distribution of the accelerated charged particles, 
even without overestimation of the photon density, gives 
us the clear indication that we should not restrict our- 
selves to the use of realistic field configurations alone. 
By not taking into account SY and/or IC losses in rel- 
ativistic particle acceleration processes in reconnection 
zones the particle pile-up at the upper energy limit is 
completely neglected and the maximum energy of the 
particles is incorrectly estimated. This shows that the 
consideration of loss processes should be obligatory for 
the determination of quantitative results in the relativis- 
tic case. In this context we briefly note that particle 
acceleration by 3D-reconnection has considerable obser- 
vational consequences. As is obvious from our simula- 
tions 3D-reconnection produces power-law energy dis- 
tributions, which would emit power-law radiation spec- 
tra, i.e. decreasing radiation flux with increasing fre- 
quency v: Flux oc j/^". For synchrotron and inverse 
Compton emission the index of the energy distribution p 
is related with the spectral index a via p = 2a + 1 In con- 
trast the simple analytical models with a constant elec- 
tric field give mononergetic energy distributions which 
give synchrotron spectra where the received radiation 
flux increases with increasing frequency up to some crit- 
ical frequency determined by the energy of the acceler- 
ated particles (Flux-ac v^/^ for an optically thin monoen- 
ergetic source) (c.fJU). 

At the present time we do not consider any interaction 
between particles and plasma waves and we neglect losses 
due to pair production. These effect should be addressed 
in future simulation studies as well as the interplay of 
second order magnetic flelds induced by accelerated par- 
ticle beams and the ponderomotive force of the excited 
plasma fluctuations. We feel that such studies will prove 
to be very helpful to clarify the role of reconnection in a 
variety of dynamical astrophysical plasma phenomena as 
different as, e.g., extragalactic jets, TeV burst and solar 
flares. 
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